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Abstract 

It is demonstrated that under common conditions a molecular solid subject to Jahn- Teller in- 
teractions supports stable Q-ball-like non-topological solitons. Such solitons represent a localized 
lump of excess electric charge in periodic motion accompanied by a time-dependent shape distortion 
of a set of adjacent molecules. The motion of the distortion can correspond to a true rotation or 
to a pseudo-rotation about the symmetric shape configuration. These solitons are stable for Jahn- 
Teller coupling strengths below a critical value; however, as the Jahn- Teller coupling approaches 
this critical value, the size of the soliton diverges signaling an incipient structural phase transition. 
The soliton phase mimics features commonly attributed to phase separation in complex solids. 
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The inventory of electronic and optical materials with technological applications is contin- 
ually expanding. One theme that is emerging in these complex materials is the importance 
of the interplay between the electronic and lattice degrees of freedom in shaping the prop- 
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erties of the materials jjj. While some kinds of electron-lattice interactions (deformation 
potential coupling, piezoelectric coupling, etc.) have been well-studied for decades, impor- 
tant qualitative differences can emerge in situations where there are electronic or vibrational 
degeneracies, a frequent occurrence in crystals with high spatial symmetry. Jahn- Teller in- 
teractions involving the interaction between degenerate (or nearly degenerate) electronic 
and vibrational states are now thought to play a substantial role in many of the recently 
discovered materials systenrs. ■rrcluding high-tenrperature superconductors Q. alkal-doped 
fuUerides j3|, and the manganites |4|. 

Here, a Jahn- Teller model is presented that supports non-topological solitons stabilized 
by a coupling of electron and lattice degrees of freedom. The analogues of such states have 
)een coniectured to have relevance to cosmological phase transitions in the early Universe 
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3], and they have been cited as potential candidates for dark-matter p]. Thus the materials 
discussed here are experimentally accessible systems where these states may be studied. The 
soliton phase in this model provides an intrinsic mechanism for the kind of heterogeneous 
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and "phase separation" found with increasing frequency in complex 



charge clumping 
materials 

The Model Consider a molecular crystal describing three electronic bands linearly cou- 
pled to a doubly degenerate distortion. The electronic bands are constructed from an orbital 
doublet, transforming as x and y and a trivial orbital singlet; the doublet of distortions trans- 
forms as X and and it is treated as a doublet of classical fields. The following continuum 
model Hamiltonian is taken 

H = H, + Hph + HjT (1) 



where 



He= J2 d^xi-li^l'^^i^m + W6\^miii^lij„,) (2a) 

m=0,±l 

Hph= d^x{\dt(f)m\'^ + K2\(f)m\'^ + K3\(f)^f + K4\(j)ni\^) (2b) 
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-9 J2 d^<<P-m^L^o + H.c.) (2c) 
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m=±l ' 
HjT 

m=±l 



where W is the energy sphtting between the orbital doublet and the singlet, and m labels 
the states in the axial angular momentum basis. Atomic units {h = rrie = 1) are used 
throughout. 

This is a continuum version of the pseudo-Jahn- Teller model with local symmetry {A © 
E) © e, and it is similar to a recently studied continuum model that was shown to contain 
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new kinds of multiorbital polarons ( "vector polarons" ) [8|] that feature spatial variations of 
the vibronic mixing. It is demonstrated here using a variational argument that the model of 
Eq. ^admits non-topological Q-ball-like solitons P,!^ for a range of the model's parameters. 

With a doubly degenerate vibrational spectrum, the model contains the following con- 
tinuous symmetry 

(pm exp{ima)(l)m (3a) 



^rn exp{-ima)ijjm (3b) 

where m = 0, ±1. This global symmetry corresponds to conservation of axial angular 
momentum of the combined electron-lattice system. The Noether charge associated with 
this global transformation is found to be given by 

Q = ^ (fx m{ll)lj)m + i(j)radt(f)-ni) (4) 

Variational Ansatz Consider a trial form for the distortion fields with harmonic time- 
dependence 

0±i(^, t) = ^(r) exp{±iujt) (5) 

corresponding to a spatially varying amplitude and a region with a local molecular distortion 
that rotates about the symmetric configuration. Additionally the amplitude of the distortion 
^(r) is taken to be a spherical pulse in space, consistent with the Q-ball field amplitude in the 
thin wall approximation . The amplitude and radius of the pulse are treated as variational 
parameters to be determined by minimizing the energy of the configuration subject to the 
constraint of fixed Noether charge. 

In the frame co-rotating with the distortion, the Hamiltonian acquires an additional term 
responsible for rotation-electronic coupling, a term well-known in the effect of A-doubling 
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Hr = -i^Y / ^i^ln^m (6) 

m 

This additional term respects the symmetry transformations of Eq. IHl and the Noether 
charge Q is also conserved in the corotating frame. 

The Bogoliubov-de Gennes equations for the Hamiltonian system H^, + Hjt + are 

-]^V'^uo- g^{u_i + ui) = euo (7a) 



-^V^Mi + {W - uo)ui - giuo = eui (7b) 



■]^V'^u^i + {W + u)u^i-giuo = eu^i (7c) 



In the case where W ^ uj and yU (the chemical potential), the system of equations reduces 

to 

^ ^^^^ 
«-i ^ TiF^^o (8b) 

1^2 V'^^W^ /o N 



Electrons in the m = channel move in an effective square well potential that is attractive. 
The strength of the potential varies inversely with W and is proportional to the square of 
the distortion. A Fermi gas experiencing this short-range attractive potential will form an 
electron-rich region in the vicinity of the potential with Friedel oscillations in the density 
that decay away with distance (see Fig. [T]). 

Within the Thomas-Fermi (TF) approximation for the kinetic energy, an analytic expres- 
sion for the total kinetic energy of the electrons in the co-rotating frame can be found. For 
W large {W ^ g^Q,fi,uj), it is sufficient to only include the kinetic energy of electrons in 
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FIG. 1: Electron charge density profile for a Jahn- Teller soliton. Local distortions create an 
attractive potential well that enhances the electron charge density. Friedel oscillations in the 
charge density set in beyond the distorted region r > tq. 

the m = channel, and only terms 0{{u!/Wy) need be retained. Dependence on uj can be 
eliminated in favor of the Noether charge Q, using the relation obtained from Eqs. EJandlHl 

Q ^ ^-0nouV. (9) 

where no is the electron density away from the distorted region and V is the volume of the 
spherically symmetric distorted region. 
Thus, 



E{V, eo) = - ^)eoV + + ^seoV + n,eoV (10) 



1, 2^?% 2 

for the case where /i ^ g^o where is the amplitude of the distortion jl2|. 

Jahn-Teller Solitons The energy in Eq. ^Jis minimized with respect to the two varia- 
tional parameters to give V"^ = (3/QU{^o) and = — '^3/2/t4 where the effective potential 
f/(^o) = 1(^2 - 2g'^nQ/W)Q + Kg^^ ^ ^^^a p ^ 3Q^Wyi6g'^no. Thus, a stable non- 
topological soliton is found for systems with negative cubic anharmonicity, a condition com- 
monly invoked in models of the thermal expansion in solids. The amplitude of displacement 

is determined only by elastic constants. This is in contrast to other Jahn- Teller- induced 
d.p.a.ve .odels Q ^ w.e. tKe an,p.,tude of *e displacement depends also on tKe 
strength of the electron-lattice coupling g. 

Furthermore, it is noted that U{$,q) vanishes at a critical coupling gc = 



^^/W{2k2K,4 — Kl)/AnQK4 (see Fig. 2). Thus the volume occupied by the soliton V diverges 
when g = gc. As a consequence, it is seen from Eq. El that to also vanishes at gc for finite 
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FIG. 2: (a) Variational energy of Jahn- Teller soliton E versus coupling constant g. The soliton 



energy behaves as \/g^ — g^/g, thus vanishing as — > Qc- The vanishing soliton energy signals the 
onset of a structural phase transition, (b) Sketch of C/(^) for several values of g. At gc, U vanishes 
at the local minimum at = I '^s 1/2^4. 



charge, and the energy of the distorted state at g^. becomes equal to the energy of the sym- 
metric (undistorted) state. Hence the system undergoes a structural phase transition for 
g > gc- It should be noted that for values of g close to gc, quantum fluctuations in the 
distortion become important and will alter the nature of this phase transition jl5| . 

While only the volume and not the shape of the soliton alters the energy of this variational 
solution, the additional energy associated with a finite-sized wall, connecting the distorted 
region of the interior to the symmetric region exterior to the soliton, contributes a surface 
tension ^. In the absence of vibrational or electronic anisotropy, the soliton is a sphere, 
the shape with minimum surface area. 

This model can be generalized in several ways. The effects of electron-electron interactions 
might be included by the addition of a Hubbard-f/ term in Eq. |23 Within mean-field theory, 
the effects of the electron repulsion work counter to the lattice-induced attractive potential 
of Eq. ISl however, for a range of U below a critical value Uc, the soliton will continue to be 
a stable solution of the model. 

Another generalization results from enlarging the symmetry group by coupling to larger 
multiplets. Consider vibronically coupling a singlet band to a degenerate triplet of bands 
with a triplet of degenerate distortions, a continuum pseudo- Jahn- Teller model with local 
symmetry (A © T) t. The triplets support the vector representation of 5*0(3). Thus the 
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vibronic interaction is invariant under the following continuous transformations: 



m' 

(lib) 

m' 

where T^m'mi-^) J = 1 Wigner "D-matrix associated with the rotation R. 

For vibronic coupling below a critical value, the 5*0(3) symmetry is unbroken and a 
stable Jahn- Teller soliton results for a region of the parameter space. The key ingredients in 
generalizing these results to other pseudo- Jahn- Teller models are adapted from Coleman P| : 
(1) the existence of an unbroken continuous symmetry, and (2) a global minimum in f/(^)/^^ 
for non- vanishing ^. For models with classical distortions, Jahn- Teller interactions (rather 
ihan pseudo) will not work even though many are invariant under continuous symmetries 
16^ . For such systems, the symmetry is broken for arbitrarily small coupling constant, thus 
violating the first criterion. 

In summary, a continuum Jahn- Teller model was introduced that supports non- 
topological (Q-ball-like) solitons that self-localize through local electron-lattice interaction. 
Such solitons are stable and are of finite size for coupling strength below a critical value. The 
size of the soliton diverges as the coupling strength approaches this critical value, signaling 
a structural phase transition in the system. It is anticipated that such states should be 
observable in two ways: (1) by measuring time-dependent fiuctuations in atomic positions 
in molecular crystals that are close to a structural phase transition, and (2) by directly 
imaging the electron charge density. A soliton phase in a molecular crystal would exhibit 
charge inhomogeneity and concomitant localized distortions, thus serving as a new kind of 
nanoscale structure relevant to a variety of complex materials. 
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